Search for: All records

Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process (r-process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We presentνh-amr, a novel neutrino-transport general relativistic magnetohydrodynamic (νGRMHD) code, which we use to perform the first 3D globalνGRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, $\dot{M}\sim 0.1-1M_{\odot }{\mathrm{s}}^{-1}$, do not produce significant neutron-rich (Y_e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, $\dot{M}\gtrsim \text{few}{M}_{\odot }{\mathrm{s}}^{-1}$(e.g., for more centrally concentrated progenitors), can unbindM_ej ≲ M_⊙of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflowY_e; however, the finalr-process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines becomer-process factories. 

Abstract Accretion disks around both stellar-mass and supermassive black holes (BHs) are likely often warped. Whenever a disk is warped, its scale height varies with azimuth. Sufficiently strong warps cause extreme compressions of the scale height, which fluid parcels “bounce” off of twice per orbit to high latitudes. We study the dynamics of strong warps using (i) the nearly analytic “ring theory” of Fairbairn & Ogilvie, which we generalize to the Kerr metric, and (ii) three-dimensional general-relativistic hydrodynamic simulations of tori (“rings”) around BHs, using theH-AMRcode. We initialize a ring with a warp and study its evolution on tens of orbital periods. The simulations agree excellently with the ring theory until the warp amplitude,ψ, reaches a critical valueψ_c. Whenψ > ψ_c, the rings enter the bouncing regime. We analytically derive (and numerically validate) that ${\psi }_{\mathrm{c}}\approx {(r/r_{\mathrm{g}})}^{-1/2}$in the non-Keplerian regime, wherer_g = GM/c²is the gravitational radius, andMis the mass of the central object. Whenever the scale height bounces, the vertical velocity becomes supersonic, leading to “nozzle shocks” as gas collides at the scale height minima. Nozzle shocks damp the warp within ≈10–20 orbits, which is not captured by the ring theory. Nozzle shock dissipation leads to inflow timescales 1–2 orders of magnitude shorter than unwarpedαdisks, which may result in rapid variability, such as in changing-look active galactic nuclei or in the soft state of X-ray binaries. We propose that steady disks with strong warps may self-regulate to have amplitudes nearψ_c. 

Abstract The observability of afterglows from binary neutron star mergers occurring within active galactic nuclei (AGN) disks is investigated. We perform 3D GRMHD simulations of a postmerger system and follow the jet launched from the compact object. We use semianalytic techniques to study the propagation of the blast wave powered by the jet through an AGN disk-like external environment, extending to distances beyond the disk scale height. The synchrotron emission produced by the jet-driven forward shock is calculated to obtain the afterglow emission. The observability of this emission at different frequencies is assessed by comparing it to the quiescent AGN emission. In the scenarios where the afterglow could temporarily outshine the AGN, we find that detection will be more feasible at higher frequencies (≳10¹⁴Hz) and the electromagnetic counterpart could manifest as a fast variability in the AGN emission, on timescales less than a day. 

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $$t_{\rm gen}$$, and one for the large-scale scale advection, $$t_{\rm adv}$$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $$\alpha$$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $$\alpha$$-dynamos. 

Abstract Disk-fed accretion onto neutron stars can power a wide range of astrophysical sources ranging from X-ray binaries, to accretion-powered millisecond pulsars, ultraluminous X-ray sources, and gamma-ray bursts. A crucial parameter controlling the gas–magnetosphere interaction is the strength of the stellar dipole. In addition, coherent X-ray pulsations in many neutron star systems indicate that the star's dipole moment is oblique relative to its rotation axis. Therefore, it is critical to systematically explore the 2D parameter space of the star's magnetic field strength and obliquity, which is what this work does, for the first time, in the framework of 3D general-relativistic magnetohydrodynamics. If the accretion disk carries its own vertical magnetic field, this introduces an additional factor: the relative polarity of the disk and stellar magnetic fields. We find that depending on the strength of the stellar dipole and the star–disk relative polarity, the neutron star's jet power can either increase or decrease with increasing obliquity. For weak dipole strength (equivalently, high accretion rate), the parallel polarity results in a positive correlation between jet power and obliquity, whereas the antiparallel orientation displays the opposite trend. For stronger dipoles, the relative-polarity effect disappears, and jet power always decreases with increasing obliquity. The influence of the relative polarity gradually disappears as obliquity increases. Highly oblique pulsars tend to have an increased magnetospheric radius, a lower mass accretion rate, and enter the propeller regime at lower magnetic moments than aligned stars. 

Abstract The spin of a newly formed black hole (BH) at the center of a massive star evolves from its natal value due to two competing processes: accretion of gas angular momentum that increases the spin and extraction of BH angular momentum by outflows that decreases the spin. Ultimately, the final, equilibrium spin is set by a balance between both processes. In order for the BH to launch relativistic jets and power aγ-ray burst (GRB), the BH magnetic field needs to be dynamically important. Thus, we consider the case of a magnetically arrested disk (MAD) driving the spin evolution of the BH. By applying the semianalytic MAD BH spin evolution model of Lowell et al. to collapsars, we show that if the BH accretes ∼20% of its initial mass, its dimensionless spin inevitably reaches small values,a≲ 0.2. For such spins, and for mass accretion rates inferred from collapsar simulations, we show that our semianalytic model reproduces the energetics of typical GRB jets,L_jet∼ 10⁵⁰erg s⁻¹. We show that our semianalytic model reproduces the nearly constant power of typical GRB jets. If the MAD onset is delayed, this allows powerful jets at the high end of the GRB luminosity distribution,L_jet∼ 10⁵²erg s⁻¹, but the final spin remains low,a≲ 0.3. These results are consistent with the low spins inferred from gravitational wave detections of binary BH mergers. In a companion paper by Gottlieb et al., we use GRB observations to constrain the natal BH spin to bea≃ 0.2. 

Abstract Black hole (BH) spin can play an important role in galaxy evolution by controlling the amount of energy and momentum ejected from near the BH into the surroundings. We focus on radiatively inefficient and geometrically thick magnetically arrested disks (MADs) that can launch strong BH-powered jets. With an appropriately chosen adiabatic index, these systems can describe either the low-luminosity or highly super-Eddington BH accretion regimes. Using a suite of 3D general relativistic magnetohydrodynamic simulations, we find that for any initial spin, an MAD rapidly spins down the BH to the equilibrium spin of 0 <a_eq≲ 0.1, very low compared toa_eq= 1 for the standard thin luminous (Novikov–Thorne) disks. This implies that rapidly accreting (super-Eddington) BHs fed by MADs tend to lose most of their rotational energy to magnetized relativistic outflows. In an MAD, a BH only needs to accrete 20% of its own mass to spin down froma= 1–0.2. We construct a semi-analytic model of BH spin evolution in MADs by taking into account the torques on the BH due to both the hydrodynamic disk and electromagnetic jet components, and find that the low value ofa_eqis due to both the jets slowing down the BH rotation and the disk losing a large fraction of its angular momentum to outflows. Our results have crucial implications for how BH spins evolve in active galaxies and other systems such as collapsars, where the BH spin-down timescale can be short enough to significantly affect the evolution of gamma-ray emitting BH-powered jets. 

Abstract The conventional accretion disk lore is that magnetized turbulence is the principal angular momentum transport process that drives accretion. However, when dynamically important large-scale magnetic fields thread an accretion disk, they can produce mass and angular momentum outflows, known as winds,that also drive accretion. Yet, the relative importance of turbulent and wind-driven angular momentum transport is still poorly understood. To probe this question, we analyze a long-duration (1.2 × 10⁵r_g/c) simulation of a rapidly rotating (a= 0.9) black hole feeding from a thick (H/r∼ 0.3), adiabatic, magnetically arrested disk (MAD), whose dynamically important magnetic field regulates mass inflow and drives both uncollimated and collimated outflows (i.e., winds and jets, respectively). By carefully disentangling the various angular momentum transport processes within the system, we demonstrate the novel result that disk winds and disk turbulence both extract roughly equal amounts of angular momentum from the disk. We find cumulative angular momentum and mass accretion outflow rates of $\dot{L}\propto r^{0.9}$and $\dot{M}\propto r^{0.4}$, respectively. This result suggests that understanding both turbulent and laminar stresses is key to understanding the evolution of systems where geometrically thick MADs can occur, such as the hard state of X-ray binaries, low-luminosity active galactic nuclei, some tidal disruption events, and possibly gamma-ray bursts. 

Abstract Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius,R_B= 10³R_g. The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outsideR_Band suppress BH accretion to 1.5% of the Bondi rate, ${\dot{M}}_{\mathrm{B}}$. Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% of ${\dot{M}}_{\mathrm{B}}$reach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2% ${\dot{M}}_{\mathrm{B}}$. Because the BAD and RAD states tangle up the jets and destroy them well insideR_B, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies. 

Abstract Accreting supermassive black holes (SMBHs) produce highly magnetized relativistic jets that tend to collimate gradually as they propagate outward. However, recent radio interferometric observations of the 3C 84 galaxy reveal a stunning, cylindrical jet already at several hundred SMBH gravitational radii,r≳ 350r_g. We explore how such extreme collimation emerges via a suite of 3D general relativistic magnetohydrodynamic simulations. We consider an SMBH surrounded by a magnetized torus immersed in a constant-density ambient medium that starts at the edge of the SMBH sphere of influence, chosen to be much larger than the SMBH gravitational radius,r_B= 10³r_g. We find that radiatively inefficient accretion flows (e.g., M87) produce winds that collimate the jets into parabolas near the black hole. After the disk winds stop collimating the jets atr≲r_B, they turn conical. Once outsider_B, the jets run into the ambient medium and form backflows that collimate the jets into cylinders some distance beyondr_B. Interestingly, for radiatively efficient accretion, as in 3C 84, the radiative cooling saps the energy out of the disk winds; at early times, they cannot efficiently collimate the jets, which skip the initial parabolic collimation stage, start out conical near the SMBH, and turn into cylinders already atr≃ 300r_g, as observed in 3C 84. Over time, the jet power remains approximately constant, whereas the mass accretion rate increases; the winds grow in strength and start to collimate the jets, which become quasi-parabolic near the base, and the transition point to a nearly cylindrical jet profile moves outward while remaining insider_B.